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We present an interpretation of tlie recent neutron data on tlie commensurate and incommensu- 
rate spin fluctuations found in YBa2Cu306+j; based on a special configuration of the electronic 
dispersion and intervention from the d^2_y2-wave superconducting phase. The observed switch over 
between the commensurate and incommensurate fluctuation spectra at the change of frequency or 
temperature is naturally accounted within this scenario. 



Recent inelastic neutron scattering (INS) experiments 
on YBa2Cu306+i, (YBCO) compounds show a lot of new 
boggling data in spin dynamics. Since magnetic fluctua- 
tions in the high- Tc cuprates are long thought to be inti- 
mately related to their mechanism of superconductivity, 
it has attracted very much attention. 

We briefly review the INS experiments which mea- 
sure spin susceptibility ln\Xs on YBCO. The INS in low- 
temperature superconducting (SC) phase optimal-doped 
YBCO shows a peak at frequency loq{T—Q) = 41meV 
and momentum ~ (tTiTt) ||^. The frequency width 
of the peak is within lOmeV and momentum width is 
within O.Itt - O.Svr [U|. In underdopcd YBCO, the res- 
onance frequency is slightly softened, and some damp- 
ing is developed both in the frequency and momentum 
space [Q . When the temperature is raised, the resonance 
frequency loo{T) is softened only very little and inten- 
sity slightly suppressed before it abruptly disappears at 
the transition to normal state j2j . Recently it was found 
that the incommensurate fluctuation, which was found 
previously only in La2-KSra;Cu04 (LSCO) compounds, 
also exist in YBCO HJ^. The spin fluctuation in low- 
temperature YBCO at frequencies either lower or higher 
than wo(0) is incommensurate in nature, characterized 
by four peaks at (tt, tt ± 6) and (tt ± 5, n) as in the case 
of LSCO. At frequencies just below ijJo(O) the deviation 
of the incommensurate peak from Q^r is also found to 
be slightly less than those at lower frequencies, where in 
contrast the deviation grows with frequency at w > cjo (0) 
There is also evidence that at frequency substan- 
tially below wo(0), raising the temperature from below 
to high above Tc brings the fluctuation from incommen- 
surate to commensurate Contradictory, in a similar 
underdoped YBCO, a weak incommensurate structure is 
found recently by Aral et Q at a similar high temper- 
ature regime. 

On the other hand, a similar commensurate peak was 
also found very recently in BiaSraCaCuzOg (BSCCO) 0. 
Though less definite, an incommensurate fluctuation also 
occurs in these compounds at lower frequencies |^ 

Since in earlier times commensurate fluctuation was 
specific in YBCO while incommensurate in LSCO, the 
fluctuations were believed to be mutually unrelated and 
hence theoretical approaches to them were independent. 



The recent found interrelation thus sets a new constraint 
on the candidate theory of these systems. 

The incommensurate peaks are believed as, due to 
Fermi surface nesting in one form or another |P,p^. 
Whereas for the commensurate peak, various possibili- 
ties were proposed. Some focused on its peak feature 
in frequency space while others gave more com- 

prehensive treatments including the overall behavior of 
luiXs in momentum space |]l0| , |l5| -|l9| . Because of the re- 
cent discovery of incommensurate fluctuation in YBCO, 
unified treatments of the fluctuations were also put for- 
ward 1 10 1^. A common feature of these theories is 
that most of them incorporate some interaction effect 
into the susceptibility via Random-Phase- Approximation 
(RPA) in some form to result in a collective mode at Q^r 
|lO|,yjl5|-0,|l|l. In this Letter we address the mentioned 
observations to a property of the band structure in con- 
junction with the intervention from the d^2_y2-wave SC 
phase. Thus in contrast to the collective mode via RPA 
scenarios, our view of the 41meV peak is based on the 
bare susceptibility. Furthermore we also include the dis- 
cussion of the switch-over between commensurate and 
incommensurate fluctuations as driven by temperature, 
which so far absent from the literature. 

Our approach was initiated by recent angle resolved 
photoemission spectroscopy (ARPES) measurement on 
YBCO and BSCCO. The measurement shows the exis- 
tence of an extended van Hove singularity (VHS) at the 
M points [(0, ±7r) and (±7r, 0)] of the Brillouin zone situ- 
ated within 20meV below the Fermi level ^j^. The 
authors relates this to the occurrence of higher Tc of 
these compounds as enhanced by the VHS |^ , compared 
to NdCeCuO compounds where the VHS is found hun- 
dreds of meV below the Fermi level |2^]. On the other 
hand, both YBCO and BSCCO were found to exhibit 
commensurate fluctuation, and incommensurate fluctua- 
tion at lower frequencies. We believe that these common 
features are not accidental. Therefore in this paper we 
start out from an electronic dispersion with Fermi level 
near the VHS at M to describe the INS experiment on 
YBCO. Previously there were also interpretations of the 
commensurate peak based on the VHS effect. But the au- 
thors considered the RPA-type susceptibility in bilayer 
systems and a s^'^-wave gap [|7j. A mono-layer 
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VHS scheme closer to ours also exist but the discussion 
on the bare susceptibility was limited to zero temperature 
and specific to commensurate peaks Q. 

Our description of the spin fluctuations in YBCO is 
in good agreement with experiments. The absence of 
commensurate magnetic response in LSCO is probably 
due to the remoteness of Fermi level from the VHS at M 
as it has a lower Tc- 

We start with the single-band bare spin susceptibility 
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where f{Ei^) is the Fermi function and = (e^ + A^)i 
is the quasiparticle excitation spectrum with Ek and Ak 
the electronic dispersion and superconducting gap re- 
spectively. We use a two-dimensional tight-binding elec- 
tronic dispersion 



Ek — — 2i(cos fcj, + cos ky) — At cos fc^ cos ky 



where t and t' are the nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) hopping respectively. We have 
fixed t' = — 0.178i throughout this paper and fi is the 
chemical potential that depends on doping. With YBCO 
in mind, the gap is taken as Ak = A(T)(cosfca; — 
cosfcy)/2. The maximum gap at zero temperature is 
taken as A(0) = 0.3t and this point will be further 
elaborated. The numerical integration width is taken as 
5 = omit and the slicing is 2000 x 2000. 

The chosen electronic dispersion is meant to phe- 
nomcnologically describe the empirical relation between 
the Fermi level and the VHS. The chosen NNN hopping 
t = — 0.178t has a Fermi surface crossing M at doping 
level a; = 0.15 [Fig. 1(b)], which is here nominally taken 
as the optimal doped case. This is to simulate a most 
intense effect of the VHS at an optimal doped compound 
as to result at a highest T^. We choose a; = 0.10 through- 
out this paper unless otherwise stated, as to show that 
the commensurate peak exist even when the Fermi sur- 
face is not exactly situated at the VHS but sufficiently 
close to it. And we believe that the only crucial point 
is the small difference in energy between the VHS and 
Fermi level but not their distance in momentum space 
as seen in the ARPES experiment. Whatever the origin 
of the VHS is also out of our concern. Furthermore we 
adopt a rigid band approximation. As far as the chosen 
dispersion has a Fermi surface out-bowed from the origin 
of the Brillouin zone as seen in the ARPES experiment, 
the main conclusion of our discussion is believed to valid. 
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FIG. 1. (a) Zero-temperature Imx2(Q7r, t^) at doping 
X = 0.05, 0.15, and 0.25. When making comparison of the res- 
onance frequencies in an absolute energy scale, the empirical 
fact of A(0) in underdoped YBCO could be smaller should be 
taken into account. The broadening of the peak upon under- 
doping is discussed in the text, (b) Fermi surfaces correspond 
to the doping levels. The Q,r transition corresponding to the 
leading edge of the peak at a; = 0.05 is shown. 



In Fig. 1(a), the frequency dependence of zero- 
temperature Imx°(Q7r, ti-") at several doping levels rang- 
ing from underdoped to overdoped are shown to exhibit 
a peak near 2A(0). Imx°(Q7r,w) is gapped out at low 
frequencies since Q.^- spans regions of maximal gap. It 
shows a high intensity at a small window around 2A(0) 
due to the enhancement by the coherence factor on tran- 
sitions between opposite gap phase regions, and is most 
intense when the frequency is just enough to open up the 
transition [which is near 2A(0)]. Furthermore, these re- 
gions are situated near M. Thus the peak arises due to 
a co-enhancement from the VHS and coherence factor. 
Note that since the gap is written in term of the NN- 
hopping, which is generally different from system to sys- 
tem, the absolute resonance frequencies and linewidths 
of different doped cases could not be directly compared. 
But if we phenomenologically input the experimental fact 
that cJo(O) is only moderately decreased at underdoping, 
the peak could be argued as always be broadened at un- 
derdoping. We give the argument as follows. Denot- 
ing F as the half-height width of the peak. At doping 
X — 0.05, with choices of A(0)/t within 0.2~0.5, the ra- 
tio r/wo(0) runs within 0.14~0.07. While at a: = 0.15, 
the same choice of A(0)/t give F/wq within 0.05~0.04 (in 
both cases a smaller A(0)/f gives a larger F/ci;o(0)). It 
is seen that if the decrement of is within 30 percent, 
the broadening in an absolute energy scale always exist 
regardless of the choice of A(0)/i in individual systems 
within a broad range. 

The broadening of the peak in frequency space is due 
to a mismatch of both enhancement effects. As the com- 
pound is underdoped, the Q^r transition on the Fermi 
surface shifts away from the full-gap regions [see Fig. 1(b)] 
and leads to a down shift of the leading edge of the peak. 
On the other hand the VHS remains to affect only when 
the frequency is sufficient to open up the transition to 
the M points. In the overdoped case, Imx°(Q^,w) is 
suppressed because the Fermi surface has shrunk to a 
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FIG. 2. The frequency evolution of Imxs in momen- 
tum space at T = and x = 0.10. (a) u = 0.3A(0); 
(b) a; = l.OA(O); (c) oj = 2.025A(0) [i.e. tJo(O)]; (d) 
liJ = 2.2A(0). (e) Scans along q = (O.Ttt, 7r)^(7r, tt) at fre- 
quencies w/A(0) = 1.0, 1.3, and 1.7. The inset shows the 
Lj-dependence of the momentum location for the tips of the 
peaks, (f) Scans along q = (O.Stt, 0.57r)^(7r, 7r)^(1.57r, tt) at 
c^/A(0) = 2.2,2.6, and 3.0. 



size that Q^r is too long to connect the Fermi surfaces, 
which all transitions are bound on [Fig. 1(b)]. 

Imx2(q, Wo) in momentum space is indeed a peak at 
with a reasonable width [Fig. 2(c)]. In our calculation 
the broadening effect of underdoping is more obvious in 
frequency space than in momentum space. Another im- 
portant point to note is that apart from the nearness of 
the Fermi level to the VHS, the existence of the peak 
also demand the existence of a d^^ -y2-w&Ye gap, but not 
a mono-layer s-wave gap. 

Fig. 2 shows the zero-temperature Imx°. It consist of 
node-to-node incommensurate peaks at quasielastic fre- 
quencies, Fermi-surface-nesting incommensurate peaks 
at intermediate frequencies, commensurate peak at fre- 
quency wo(0) near 2A(0), and incommensurate peaks at 
frequencies over ijJq{0) [Fig.2(a)-(d)]. The origin of the 
node-to-node and nesting incommensurate peaks are well 
accounted in the literature [p[p^ . In our scenario such se- 
quential emergence of the peaks naturally occurs in con- 
trast to the RPA scenarios where the resonance frequency 
of the commensurate peak could be sensitive to the choice 
of the interaction strength. 

There is a decrease of the incommensurateness of the 
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FIG. 3. (a) Temperature dependence of lm'x'l{Q,T^,ijj) at 
uj = 2.025A(0) [i.e. t^o(O)], 1.5A(0), and l.OA(O). (b) 
c^-dependence of Imx° at and T = 0,0.77rc, 0.92rc, and 
Tc. The doping is fixed ai x — 0.10 and the temperature 
dependence of gap A(r) is shown in the inset of (a). 



nesting incommensurate peaks at frequencies just below 
cjo(O) [Fig. 2(e)]. This is due to a dynamical nesting ef- 
fect. Since the Fermi surface is out-bowed, the opened 
up energy contour E'k = uj/2 which determines the possi- 
ble transitions [|I^ has a smaller (larger) curvature above 
(below) the Fermi surface. Thus local nesting is better 
between contour above the Fermi surface and this deter- 
mines the location of the tip of the peak. As w increases, 
the contour surges away from the Fermi surface and that 
pushes the nesting peaks towards Qtt- 

At frequencies over a;o(0), due to the filling up of the 
valley at Q^r a hardly discernible incommensurate feature 
is seen [Figs. 2(d) and 2(f)]. It is characterized by incom- 
mensurate peaks displaced from Qtt along both the zone 
diagonals and edges. The incommensurateness grows 
with frequency in this regime Fig. 2(f)] and that agrees 
qualitatively with experiment [^. Further comment will 
be made on this point at the end of this Letter. 

In our following discussion of the temperature evolu- 
tion, we take a convenient empirical relation between the 
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gap and temperature A(r)/A(0) = [1 - (T/Tc) 
inset in Fig. 3(a)] and a typical ratio 2A(0)/Tc = 8. 

In Fig. 3(a), Imx° at (Q7r,^o(0)) is seen to be sup- 
pressed by temperature as in experiments [^. At tem- 
perature below Tc, the main cause of the suppression is 
the shift away of the mode frequency luo{T) from luo{0), 
the diminishing of the peak intensity is actually less rapid 
[Fig. 3(b)]. For temperatures T < 0.8Tc, the mode fre- 
quency of the commensurate peak is softened a little from 
luq (0) before abruptly softened to zero and disappears at 
Tc. The softening is merely a manifestation of the dimin- 
ishing gap and owing to Tc which is a small energy scale 
compared to the bandwidth, the peak feature is main- 
tained well at T < Tc. 

While the temperature evolution of the commensu- 
rate peak at a; = ujo{0) has been frequently discussed 
|12| , |l3 16 1, the temperature evolution of the incommen- 
surate peaks at w < ti;o(0) has been overlook. Neverthe- 
less recent experiments show that this is not at all trivial 
[^,^. At low frequency and temperature, Imx^ at Q^r 
is gapped out and leads to a clear appearance of the in- 



3 




FIG. 4. Finite-temperature Im^s in momentum space at 
an intermediate frequency uj — l.OA(O). A(0) = 0.3t and 
2A(0)/rc = 8 throughout, (a) T = 0.955Tc and x = 0.10; (b) 
T = Tc and a; = 0.05. 



ness of the Fermi level to the VHS at M points and 
coherence effect from the dx2-y2-weLve SC phase. The 
results describe naturally the frequency and temperature 
dependence of the commensurate peak, the switch over 
between the incommensurate and commensurate spectra 
at the change of frequency or temperature, and the de- 
crease [increase] of the incommensurateness at frequency 
u! < 2A(0) [u! > 2A(0)]. Furthermore, the broadening in 
frequency space of the commensurate peak at underdop- 
ing is also accounted semi-phenomenologically. 

We thank JH Kok and YT Huang for helpful conversa- 
tions and acknowledge the support from NSC of Taiwan 
under grant No.89-2112-M-003-009. 



commensurate peaks. As the temperature is raised, the 
particular frequency could equal to the temperature de- 
pendent resonance frequency and the fluctuation is then 
converted to commensurate at a cj < ujq(0) and T 
[compare Fig. 2(b) and Fig. 4(a)]. Plotting Imx° at the w 
and Q^n- against the temperature shows the emergence of 
the peak at a temperature just below Tc [Fig. 3(a)]. As 
the temperature is raised further into the normal state, 
the fluctuation reverts to a very weak incommensurate 
feature ||2^ . Experimentally it may well escape from de- 
tection and seen as a broad commensurate structure. At 
further undcrdoping, say x — 0.05, even a well distin- 
guished commensurate feature may appear in the normal 
state [Fig. 4(b)] and it is suppressed as the temperature 
increases. Note that the intensity of this high tempera- 
ture commensurate structure is of the same order of mag- 
nitude as its low temperature incommensurate counter- 
part as seen in the experiment In our case, the valley 
at appears and leads to a pronounced incommensu- 
rate structure in the normal state only in the overdoped 
case [ p3| . This may resolve the apparent contradiction 
between the observation of a broad commensurate struc- 
ture ||^ , and the observation of a weak incommensurate 
structure [^ at similar temperature (> T^) and frequency 
(< loq) regimes in an underdoped YBCO. In our calcula- 
tion these weak incommensurate structures are deviated 
from along the zone edges. 

We make a last comment on our result. If a Hubbard 
repulsion of l.Qt is incorporated via a RPA spin suscepti- 
bility, structure at the vicinity of Q.^- will be more promi- 
nent than structures away from Q^r. Furthermore, at 
zero temperature a very sharp and intense resonance at 
Qtt occurs at ujo — 1.8A(0). The incommensurate struc- 
ture at w > cjo is prominently away from Q.^^ along the 
Brillouin zone edges. If a smaller repulsion Q.bt is taken, 
every structure resembles those of the bare susceptibil- 
ity. Since we have no reason to believe that the repulsion 
should be finely 1.0^, we believe that such a tuning is a 
less probable description of reality. 

In conclusion, we have provided an unified explanation 
of the INS experiments on YBCO based on the close- 
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